Optimization of a Langmuir-Taylor detector for lithium 
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This paper describes the construction and optimization of a Langmuir-Taylor detector for lithium, 
using a rhenium ribbon. The absolute detection probability of this very sensitive detector is mea- 
sured and the dependence of this probability with oxygen pressure and surface temperature is 
studied. Sources of background signal and their minimization are also discussed in details. And 
a comparison between our data concerning the response time of the detector and literature values 
is given. A theoretical analysis has been made: this analysis supports the validity of the Saha- 
Langmuir law to relate the ionization probability to the work function. Finally, the rapid variations 
of the work function with oxygen pressure and temperature are explained by a chemical equilibrium 
model. 



I. INTRODUCTION 

Very sensitive detectors for neutral atoms, which were 
first necessary for Rabi experiments (see reviews in ref. 
jl], 0|), still govern the feasibility of many experiments. 
Two different detectors fulfill the requirements of high 
efficiency and low background: 

• the Langmuir-Taylor detector & |], || is based on 
the surface ionization process. For ground state 
atoms, the ionization probability is large when the 
ionization potential is low. Therefore, this detec- 
tor is mostly used with alkali atoms. For thermal 
atoms, the detection probability is almost indepen- 
dent of their velocities. The detection of a few 
atoms per second is feasible, if the produced ions 
are detected with an electron multiplier. 

• the laser induced fluorescence detector, which is 
commonly used in cold atom experiments, has a 
detection probability close to 100% for slow atoms. 
However, its detection probability decreases rapidly 
with the atom velocity and a large efficiency is very 
difficult to achieve for thermal atoms. Finally, laser 
stray light is usually the dominant source of back- 
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ground and limits the detection of a very low atomic 
flux. 

Therefore, the Langmuir-Taylor detector is probably 
the best detector for thermal alkali atomic beams. Op- 
timization of this detector is more difficult in the case of 
lithium, because lithium has the highest ionization po- 
tential of the alkali atoms. In the present paper, we give 
a detailed description of the detector in this case. More 
precisely, after the introduction, we describe the detector 
principle and its design in part 2. In part 3, we present 
our measurements of the detection probability, which de- 
pends on the degree of oxidation of rhenium surface Q . 
The various sources of background signal and their min- 
imization are described in part 4, while, in part 5, we 
discuss the detector response time. In part 6, these ob- 
servations are rationalized by a theoretical modeling of 
the surface ionization process. Two appendices present 
complementary information: appendix A discusses the 
relation between heating current and rhenium ribbon 
temperature and appendix B briefly analyzes literature 
data concerning the dependence of rhenium work func- 
tion with surface oxidation |?], 0, [h], |llj . 

All the information needed to optimize this detector 
is thus collected and analyzed in the present paper. A 
recent paper by F. Stienkemeier et al. |l2| has described 
the Langmuir-Taylor detector, using a rhenium surface, 
applied to detect various atoms (including lithium) at- 
tached to helium droplets. Our measurements and anal- 
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ysis are largely complementary to those of this work. 



II. DETECTOR PRINCIPLE AND DESIGN 

A. Detector principle 

Surface ionization of an atom A occurs if the ioniza- 
tion potential I of the atom is comparable to the work 
function $ of the metal. The atom is then emitted as 
a positive ion A + with a probability P + and as a neu- 
tral atom with the probability (1 — P+). The probability 
P+ is usually assumed to be given by the Saha-Langmuir 
law: 



l + ^e X p(^|) 

where go and g+ are the statistical weights of the ion 
and atom ground states (in the case of alkali atoms, 
go/ 9+ = 2). The validity of this law is discussed in part 
6. The wire temperature must be high, typically 1500 
K or larger, not because of its influence on the ionized 
fraction P + , but to reduce the ion residence time on the 
surface. This residence time r is given by: 

r = t exp (E ads /k B T) (2) 

Here E a ds is the ion adsorption energy (typically a few 
eV) and To should be close to the vibrational period of the 
ion near the surface (near 10 -13 s, see table 2). Finally, 
the ionization probability is independent of the initial 
kinetic energy as long as this energy is smaller than or 
comparable to the adsorption energy, because the resi- 
dence time on the surface is sufficient to insure thermal 
equilibrium with the surface. 

Because of the high ionization potential of lithium, 
In — 5.392 eV, we need a metal with a large work func- 
tion Among refractory metals, the highest values of 
the work function are for platinum ($ = 5.77 eV), irid- 
ium ($ = 5.70 eV) @, rhenium ($ = 4.96 eV) and 
tungsten ($ = 4.54 eV) JTlj ] . Moreover, oxidation in- 
creases the tungsten work function ($ ss 6 eV following 
N. F. Ramsey @, $ 5.9 - 6.2 eV following H. Pauly 
and J. P. Toennic s [Til ) and a similar effect occurs for 
rhenium || 0, ||[ |ll|. In appendix B, we discuss the 
dependence of rhenium work function with temperature 
and oxygen partial pressure: this information supports 
strongly the idea that rhenium is probably the best ma- 
terial for a lithium detector. The work function is also 
a function of the alkali surface coverage |7|, ||, [Tj| . This 
effect is very important if one uses a high flux especially 
when the residence time is long (i.e. at low temperature) 
but it is usually negligible for a high sensitivity detector 
dedicated to very low atomic flux. 
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FIG. 1: Schematic drawing of the detector design: 1 repeller 
plate; 2 endview of rhenium ribbon; 3 electrostatic lens; 4 
channeltron. The arrow gives the scale and a typical choice 
of potentials is also indicated. 



B. Detector design 

The central part of the detector is a rhenium ribbon. 
The ribbon thickness a must be quite small, to facili- 
tate resistive heating and cleaning (see part 4), while the 
width b is fixed by experimental needs. For example, we 
use commercially available ribbons with a = 0.0008" = 
20.3 fxm and b = 0.040" = 1016 fim provided by A. D. 
Mackay or a — 30 /im and b — 760 fj,m provided by 
Goodfellow. The ribbon is a few centimeters long. It is 
heated by circulating a DC current of a few amperes, nec- 
essary to reach the operating temperature T ~ 1500 K. 
In appendix A, we discuss the relation between the cur- 
rent and the ribbon temperature. Higher temperatures 
(T > 2200 K) are needed during the cleaning process. At 
T = 2000 K, the power per unit length is about 5 W/cm 
and this heat load may induce some outgassing of the 
vacuum tank. 

The ribbon expansion |l^, |l6| when going from ordi- 
nary temperature to a temperature near 2000 K is close 
to 1%. Usually, the ribbon is kept straight by a soft 
spring, with an applied force less than 1 Newton (for 
data concerning the mechanical strength of rhenium, see 
ref. jl^, [l6|]). In our arrangement, the rhenium ribbon is 
spot-welded at both ends to tantalum sheets. 

The ions emitted by the wire are collected and focussed 
on the entrance funnel of a channeltron, by a simple ion 
optics, as shown in figure ^. The ion trajectories have 
been calculated with the SIMION software (ItJ and this 
calculation serves to define the various electrical poten- 
tials, which are further optimized by maximizing the ion 
signal. The repeller plate potential must not exceed a 
few volts, otherwise the electrons emitted by the hot wire 
may ionize the background gas, thus contributing to the 
background signal. A fast counting electronics is used to 
convert the channeltron pulses into a signal expressed in 



3 



counts per second. 

III. DETECTION PROBABILITY 

The principle of our measurement is to use an effusive 
atomic beam of lithium. As the theory of such a beam 
is perfectly understood we can calculate the atomic 
flux dN/dt reaching the ribbon: 



where the beam intensity X is given by T = 
nLiVmeanO-* I {^) (p-Li is the lithium density in the oven, 
= ^/8fcsT firm is the mean velocity inside the 
oven, a* is the area of the oven exit hole) and Af2 is 
the solid angle of the rhenium wire seen from the oven. 
In some experiments, we used a piezoelectric slit, whose 
width is tunable under vacuum, to reduce this solid an- 
gle, thus verifying the detector linearity with the atomic 
flux. Under these conditions, the signal remains in the 
linearity domain of the channeltron as long as the lithium 
oven temperature does not exceed 673 K. 

Several equations relating the lithium vapor pressure 
to the temperature appear in the literature [|l8| [H| . Nes- 
meyanov [jl9| considers that the most reliable data cover- 
ing the 735—915 K range are represented by the equation: 

logw (pu) = 8.012 - 8172/T (4) 

with the pressure in Torr. We use the perfect gas approx- 
imation to convert pressure to density by riLi = p^/hsT. 

As we must extrapolate this equation, we made a laser 
absorption experiment with the same atomic beam to 
test this extrapolation. Because of laser saturation and 
optical pumping, the theoretical description of laser ab- 
sorption is not very simple and we use it here just as a 
relative measurement. The atomic density n in the beam 
is proportional to the absorption: 

n cx ln(/ t // ) « (/ - h)/h 

where Iq and It are the incident and transmitted laser 
intensities. This experiment was made in the tempera- 
ture range 673 — 777 K. To calibrate this measurement, 
we use the vapor pressure at the highest temperature of 
our study. The results of the absorption experiment are 
shown in figure ^ which plots the beam intensity X thus 
deduced. These results prove that the lithium pressure 
law can be safely extrapolated, for T > 673 K. 

On the same figure, we have plotted the beam intensity 
I deduced from our measurements of the ion signal, as if 
the detection probability was equal to 100%. These re- 
sults follow the same slope as the pressure law in this log- 
arithmic plot, further supporting the extrapolation of the 
pressure law. During this experimental run, the rhenium 
ribbon temperature was T — 1590 K and the residual 
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FIG. 2: Beam intensity X as a function of the oven tem- 
perature T. The plot presents the logarithm of Xy/T versus 
1/T (T in Kelvin) so that the theoretical curve deduced from 
equation Q is a straight line. The triangles represent our 
absorption measurements, calibrated at the highest temper- 
ature using equation (^). The squares represent the same 
quantity deduced from the Langmuir- Taylor detector signal, 
assuming a 100% efficiency. The absorption data prove that 
extrapolation of the pressure law is reasonably accurate. The 
Langmuir- Taylor detector signals also support the validity of 
the extrapolation of this pressure law. Finally, we deduce 
from this plot the detection probability 32 ± 5 %. 



gas pressure in the detector chamber was 3 x 10 -8 mbar 
(corresponding to po 2 ~ 6 x 10 -9 mbar if the residual gas 
is mostly air). From this plot, we get the detection prob- 
ability D + = 32 ± 5%. This error bar does not include 
the remaining uncertainty on the lithium vapor pressure. 

In a second series of experiments, we measured the 
detection probability as a function of surface temper- 
ature and oxygen pressure: the flux of lithium atoms 
impinging on the rhenium ribbon was estimated to be 
5.75 x 10 5 atoms per second. Before each measurement, 
the wire was flashed at T = 2250 K for 2 minutes, then 
the heating current was adjusted to reach the working 
temperature and oxygen was admitted. After a stabi- 
lization period, the ion signal was measured with the 
lithium beam on and off and the background was sub- 
stracted. We have varied the rhenium ribbon tempera- 
ture in the range 1590 — 1880 K. Two measurements were 
made by introducing oxygen with a leak valve, with oxy- 
gen pressures equal to 2.1 x 10~ 7 and 5.0 x 10 -7 mbar 
(uncorrected ion gauge readings) and one measurement 
with our best vacuum (3 x 10~ 8 mbar) corresponding to 
6 x 10~ 9 mbar of oxygen, if the residual gas is mostly 
air. We thus deduce the detection probability D + , which 
is plotted in figure 0. These new values of the detection 
probability and the first measurement (D + = 32 ± 5% at 
T = 1590 K in our best vacuum) are not in very good 
agreement. This discrepancy, which illustrates the dif- 
ficulty of such absolute measurements, may be due to a 
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FIG. 3: Measured detection probability as a function of rib- 
bon temperature for three different oxygen pressures: dots for 
P02 = 5.0 x 1CP 7 mbar, squares for po2 = 2.1 x 10 -7 mbar, 
triangles for our best vacuum (3 x 10~ 8 mbar) corresponding 
to an estimated oxygen pressure P02 = 6 x 1CP 9 mbar. From 
the saturation of these curves, we estimate that the true ion- 
ization probability is obtained by multiplying the detection 
probability by 1/E+ = 1.41 (see text). We have also repre- 
sented by a cross the value of D+ = 32% at T — 1590 K in 
our best vacuum deduced from the data of figure 0. 



variation of the composition of the residual gas in the de- 
tection chamber (the experiments corresponding to figure 
U have been made one year before the experiments with 
a variable oxygen pressure) . 

D + is the product of the ionization probability P+ 
by an ion counting efficiency E+. This efficiency com- 
bines several experimental factors (ion collection effi- 
ciency, channeltron efficiency including electronic thresh- 
old effects) and all the errors in the estimation of the 
atomic flux. We have no way of measuring directly E+, 
but the saturation of the detection probability with a 
maximum observed value D + = 0.66 must correspond to 
an ionization probability close to 100%, following litera- 
ture data [g, 0, g, |U|. From this remark, we get a 
good estimate of the ion counting efficiency E + = 0.71, 
this precise value being such that the corresponding P + 
values depend smoothly on the rhenium surface oxidized 
fraction / defined in Appendix B. This data is analyzed 
in part 6. 

The experiments of Stienkemeier et al. [T^ | involve 
lithium atoms attached to helium droplets and the rhe- 
nium ribbon is kept in a higher vacuum (residual pressure 
5 x 10~ 9 mbar) than in our experiment. The efficiency 
was estimated indirectly by comparison with other alkali 
atoms, with a peak value of 10 % near T — 1350 K. This 
value, substantially lower than our own, corresponds to 
a lower degree of oxidation of the rhenium surface, in 
agreement with the use of a higher vacuum. 



IV. OPTIMIZATION OF THE BACKGROUND 
SIGNAL 

The background signal has three main origins: 

• the components of the background gas with ioniza- 
tion potentials below 9 eV can be ionized on the 
rhenium surface 

• a fresh rhenium wire contains a few ppm of alkali 
atoms. When the wire is hot, these atoms diffuse 
inside the wire and reach the surface where they 
are emitted as ions 

• the background depends strongly on the oxidation 
of the rhenium ribbon and this can be explained by 
the emission of rhenium oxide ions. 

We are going to discuss now these three sources of back- 
ground signals. All these contributions to the back- 
ground signal can be emitted from any point of the rhe- 
nium ribbon surface, but the corresponding ions are de- 
tected only if they are focused on the channeltron en- 
trance funnel. A trick to reduce the background signal 
is to limit ion collection to a small part of the rhenium 
ribbon. In our case, the ion optics collects the ions emit- 
ted by the 5 mm long central part of the ribbon and this 
length cannot be much reduced, at least for our applica- 
tions. 



A. Ionization of the background gas 

Surface ionization of many organic molecules is pos- 
sible on rhenium surfaces and this technique is well es- 
tablished, as reviewed by Zandberg The ionization 
probability of any component of the residual gas de- 
pends rapidly on its ionization potential /, through Salia- 
Langmuir law ([j]). For a species with a partial pressure 
close to 10~ 8 mbar, about 5 x 10 11 molecules impinge on 
a 1 centimeter length of hot wire per second. Assuming 
a wire temperature T ~ 1500JT and an oxidized rhenium 
surface with $ « 6 eV, the corresponding contribution 
to the background remains below 10 2 ions per second if 
I — $ > 3 eV i.e. if I > 9 eV. Various types of or- 
ganic molecules have their first ionization potential near 
this value and can be ionized: a mass spectrum produced 
by ionization of the residual gas (produced by oil diffu- 
sion pumps) on an oxidized tungsten ribbon is shown in 
ref. |l4| ] and this spectrum is very dense. It is therefore 
very important to operate in a high and clean vacuum, 
because an oil-free vacuum contains mostly species with 
high ionization potentials (H2, H2O, CO...). Practically, 
the detector should be in a UHV chamber pumped by 
a turbo molecular pump or an ion getter pump. A cold 
trap at liquid nitrogen temperature may be useful to re- 
duce the vapor pressure of condensable species, but its 
use may not be necessary. 
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B. Cleaning the alkali content of the rhenium 
ribbon 



Following Goodfcllow JT5|], rhenium ribbons contain 4 
ppm of potassium in mass corresponding to 3 x 10 14 
potassium atoms per centimeter of ribbon length. At 
high temperatures, these potassium atoms diffuse inside 
the wire, reach its surface where they are emitted as ions. 
Using the diffusion equation, we get the value of this ion 
current, assuming an homogeneous initial potassium den- 
sity and a wide ribbon (a <C b): 



_ 8Nq 

7T 2 Tl 



E 



E 

p=Odd 



(5) 



where N is the total number of potassium atoms in the 
ribbon and q is the modulus of the electron charge. This 
sum of exponential is characteristic of a diffusion pro- 
cess, the diffusion mode of order p decaying with a time 
constant t p = Ti/p 2 . For long times, the decay becomes 
purely exponential, with the time constant t± related to 
the diffusion constant D by: 



7T 2 D 



(6) 
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FIG. 4: Ion current in picoamperes (logarithmic scale) as a 
function of time measured in hours (linear scale) during the 
baking process at a temperature T — 2250 K. The crosses 
represent the measurements of the total ion current, while 
the curve is the best fit to this data using equation (^). 

During the preparation of a ribbon made at a tempera- 
ture T = 2250 K, the repellcr plate (see figure 1) served to 
collect about 38 % of the emitted ions. The correspond- 
ing total ion current is plotted as a function of time in 
figure m These observations are well explained by equa- 
tion (ph. From the longest time constant ri ~ 3.0 x 10 5 
s, we deduce the diffusion constant for potassium inside 



rhenium D = 0.3 x 10" 15 m 2 s -1 at T = 2250 K. As this 
time constant is long and as it scales like a 2 , it is very 
important to choose a ribbon with a small thickness a, 
to minimize the duration of the cleaning process. 

After 6.5 hours of baking, we have briefly varied the 
wire temperature T and recorded the ion current as a 
function of temperature. This current is proportional 
to the diffusion constant, thus giving the temperature 
dependence of the diffusion constant D, usually described 
by an Arrhcnius law: 



D = D cxp(-E dlff /k B T) 

From our measurements, we get the activation energy of 
this diffusion process, Ediff — 5.36 ± 0.23 eV, and the 
prefactor value, Dq = 0.3 x 10" 



3 m 2 s 



C. Influence of rhenium oxidation on background 
signal 

While measuring the detection probability as a func- 
tion of surface temperature and oxygen pressure, we have 
recorded the background signal. We discuss here only the 
data corresponding for T = 1508 K, because the tem- 
perature dependence is not very simple. With our best 
vacuum (about 3 x 10~ 8 mbar), the background signal 
was I — 3 x 10 3 counts per second and this signal in- 
creased rapidly with oxygen pressure, reaching 1.6 x 10 4 
and 2.6 x 10 4 counts per second when the oxygen pressure 
was 2.1 x 10 -7 or 5.0 x 10~ 7 mbar respectively. Moreover, 
the noise of this background signal is substantially larger 
than Poisson noise (up to 5 times larger for 1 second 
counting periods). These two effects of a strong oxida- 
tion (rapid increase of the background and excess noise) 
decrease considerably the performance of the detector for 
a small atomic flux, even if oxidation gives an important 
gain on detection probability. 



D. Operating conditions and ultimate 
performances 

We have already explained the need for a high and 
clean vacuum. If one operates near T — 1500 K, a small 
oxygen residual pressure of the order of a few 10 -9 mbar 
is sufficient to reach a detection probability of the or- 
der of 10 — 30 % and to keep the background at a very 
low level. If the initial baking time has been sufficiently 
long, the potassium content of the rhenium ribbon may 
be sufficiently reduced to give a small background signal, 
as long as the working temperature is not too high. In 
Bielefeld, where a very large set of experiments have been 
done p^, an d because the background pressure, pro- 
duced by an ion getter pump, is below 1 x 10 -8 mbar, the 
typical background count rate is 100 counts per second. 
In Toulouse, where we have less experience and a less 
good vacuum (near 3 x 10~ 8 mbar produced by a turbo 
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pump), the typical background count rate is 1 — 3 x 10 3 
counts per second. 

Finally, in a high vacuum, surface contamination of 
the rhenium ribbon by cracking of organic molecules is 
very slow, but it is nevertheless necessary to clean the 
ribbon surface from time to time by flashing it to a high 
temperature (typically 2250 K) for a few minutes, to get 
a reproducible operation. 



V. TIME RESPONSE OF THE DETECTOR: 
RESIDENCE TIME OF LITHIUM ON A 
RHENIUM SURFACE 

The time response of a detector is very important for 
many applications. The dominant contribution for the 
present detector is the ion residence time on the surface, 
usually in the microsecond to millisecond range. The 
simplest measurement technique uses a chopped atomic 
beam Jl^]. Another method, based on the autocorrela- 
tion function of the ion current |22| , has also been be used 
and this method extends the measurement range down 
to the microsecond range, a value difficult to achieve by 
chopping atomic beams. All the data sets present in 
the literature have been fitted by equation (|^) namely 
r = r exp (E ada /k B T). 



temperature (K) 


7~exp 


r (ps) [12] 


T (us) [22] 


1525 


215 


114 


82 


1600 


75 


35 


35 



TABLE I: Ion residence time as a function of the wire tem- 
perature. Column 1 gives the temperature deduced from the 
heating current (see appendix A), column 2 our measurement 
of the residence time, columns 3 and 4 the calculated val- 
ues of r deduced from ref. jl^| and |E| respectively. The 
data collected by these authors cover the temperature range 
1200-1500 K (ref. Q) and 1600-2000 K (ref. Q) so that 
we have to extrapolate slightly their fitted laws. 

Using a chopped atomic beam, we collected two mea- 
surements, extending slightly the ranges covered by these 
works ^2| . Table 1 presents a comparison of our mea- 
surements to the values deduced from fitted laws given 
in these two papers. The agreement between these re- 
sults is rather poor and this is easily understood for the 
following reasons: 

• temperature measurement may be affected by sys- 
tematic errors (see appendix A) and any error on 
the temperature has a large effect, because the res- 
idence time t varies very rapidly with temperature 

• in the simplest theoretical model (for example see 
ref. iQ), the adsorption energy increases with the 
work function $ and the work function depends 
on the degree of oxidation of the rhenium surface. 
As different experiments test rhenium surfaces with 



different degrees of oxidation, the different values of 
the residence time for the same temperature may 
be a real effect. 

Table 2 collects the information concerning the resi- 
dence time of all the alkali on rhenium, so as to illustrate 
the trends followed by To and E a d s through this chemi- 
cal family. t is expected to be close to the vibrational 
period of the ion near the surface and this property, well 
verified in several cases, is not verified by the data con- 
cerning lithium. The very long extrapolation surely ex- 
plains the scatter on To with a correlated modification of 
the adsorption energy (as an example, the very different 
parameters of ref. and for lithium lead to the 
same value the residence time at T = 1600 K). The ad- 
sorption energy may also vary with surface oxidation, as 
discussed above, but this effect cannot explain the very 
large differences appearing in the case of lithium . 



reference 


alkali 


to in 10~ 13 s 


Eads in eV 


1 


23 


1 


Cs 


1.9 ±0.9 


2.01 ±0.04 


1 


24 


1 


Rb 


0.8 ±0.3 


2.28 ±0.03 


1 


24 


1 


K 


1.0 ±0.3 


2.33 ±0.03 


I 


12 


1 


K 


0.01 ±0.06 


2.64 ± 0.06 


[ 


24 


1 


Na 


0.2 ±0.1 


2.75 ±0.03 


1 


22 


1 


Na 


31. 


1.92 


1 


12 


1 


Na 


0.05 ± 0.04 


2.95 ±0.07 


1 


22 


1 


Li 


12. 


2.37 


1 


12 


1 


Li 


0.009 ± 0.006 


3.36 ±0.05 



TABLE II: Values of the preexponential factor to and of the 
adsorption energies E a <is of the alkalis on rhenium. Column 1 
gives the reference, column 2 the alkali, columns 3 and 4 the 
values of tq and E a d s respectively. 



VI. MODELING OF THE SURFACE 
IONIZATION PROCESS 

Surface ionization of a species A is usually described 
using the Saha-Langmuir equation (|l]). This implicitly 
assumes a thermal equilibrium between the desorbing 
particle and the surface and is thus independent of the 
characteristics of the electron transfer process between A 
and the surface. As will be shown below, it reproduces 
rather well the experimental observations. However, one 
can wonder about the corresponding microscopic view of 
the process and possible dynamical effects. The charge 
transfer between an atomic projectile and a free-electron 
metal surface is rather well understood p6[ and quan- 
titative descriptions are available (see e.g. [^?J for the 
alkali case). We can use this microscopic description to 
predict the surface ionization efficiency and compare it 
to the Saha-Langmuir prediction. 

As an alkali atom approaches a metal surface, its elec- 
tronic levels couple with the continuum of metal states, 
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resulting in a finite lifetime of the atomic levels. The 
corresponding width T gives the charge transfer rate be- 
tween the atom and the surface. At T = K, the electron 
is transferred from the atom to the surface if the atomic 
level is above the surface Fermi level and in the opposite 
direction if the atomic level is below the Fermi level. For 
a finite temperature T, electron transfer occurs in both 
directions according to the fractional population of the 
metallic states at the atomic level energy position [ p6[ . 
For a free-electron metal, the width of the alkali level 
varies approximately exponentially with the distance z 
to the surface 27 . Qualitatively, when an alkali leaves 



orption. The different contributions are summed with a 
thermal weighting factor to yield the surface ionization 
probability P+: 



the surface, the very large width at small z allows the al- 
kali charge state to reach thermal equilibrium; however, 
this is not true at large z and there exists a distance, 
called freezing distance, zp, beyond which the charge 
state of the desorbing alkali decouples from the surface 
p8| . Since the alkali level energy is a function E a (z) 
(for the ionic state, it roughly follows an image potential 
variation p7|), the charge state equilibrium value changes 
with z and the asymptotic charge state of the desorbed 
particle is different from its value on the surface and it a 
priori depends on the desorption velocity. This 'freezing 
distance' discussion provides a qualitative picture of the 
charge transfer between a projectile and a metal surface. 
Quantitatively, if we assume that the desorbing parti- 
cle follows a classical trajectory z(t), the evolution of its 
charge state is governed by a rate equation M : 



dn + 



-2Tf(E a (t),T)n+(t) 
-T(l-n + (*))(!- f(E a (t),T)) 



(7) 



where n + (t) is the ion charge fraction, f(E a (t),T) the 
Fermi function and E a (t) designates E a (z(t)). 

Surface ionization is modeled by solving equation ([?]) 
for a set of different classical trajectories, representing 
the different possibilities for a desorbing alkali. The en- 
ergies and widths in (|t]) are taken from the parameter- 
free description of the alkali free-electron metal surface 
study of ref. p7| . The initial state (z = Zi n i) close to the 
surface is taken ionic (n + = 1) and equation ^ yields 
the survival probability of the ion at a large distance for 
each trajectory and we then have to sum the contribu- 
tions from the different trajectories. In practice, we solve 
equation (Q) up to a large z distance for a set of total en- 
ergies, E, of the desorbing particle in the ionic channel. 
The z(t) trajectory introduced in equation (0) is com- 
mon to the ionic and neutral desorption channels, which 
is not appropriate for large z. To circumvent this prob- 
lem in the summation over the heavy particle energies, 
E, we assume that the charge state stabilizes around the 
freezing distance zp where we compare the local kinetic 
energy of the desorbing particle to the energy required 
for desorption in the ionic or neutral channel in order to 
decide whether desorption is energetically allowed or not. 
Ionic and neutral potential energy curves being different, 
this leads to a lower desorption threshold for neutral des- 



r Q °° P+(E) exp j-E/k B T) dE 

1 - P + " f™ h (1 - P+(E)) exp (~E/k B T) dE 



(8) 



where E is measured with respect to the ionic threshold 
and E t h is equal to E t h — E a (oo) — E a (zp). The freezing 
distance, obtained from the level width calculated in ref. 
p7j , is equal to 9.6 a for lithium; it does not vary much 
in the energy range concerned in surface ionization. 

We thus get the surface ionization probability for a 
given surface work function and temperature. Such an 
approach is indeed valid for a free-electron metal and 
should also be meaningful in the case of low adsorbate 
coverage of a metal. In the case of a low adsorbate cover- 
age on the surface, the perturbation of the charge transfer 
is usually described in terms of local and non-local effects 
(see a review on the adsorbate effects in [pl| ). The local 
effect is due to the local potential around the adsorbate 
which can strongly perturb the electron transfer in a cer- 
tain region surrounding the adsorbate. The non-local 
effect comes from the surface work function change in- 
duced by the adsorbate. An approach like the present 
one only considers the non-local aspects. Local aspects 
are mostly visible in scattering experiments which select 
specific trajectories, thus probing specific areas on the 
surface. Surface ionization a priori concerns the entire 
surface and thus should average over the local effects. For 
very large adsorbate coverages, the electronic structure 
of the surface is modified, possibly leading to an insu- 
lator layer on the surface, on which the charge transfer 
properties are different (see e.g. a review in p{|). 

We are going to describe surface ionization of lithium 
on a rhenium surface as a function of the temperature T 
and of the residual pressure, within the above approach 
and we will use the rhenium work function extracted by 
Kawano et al. jnj from their electron emission experi- 
ment (either directly the published data or the modeling 
of these results presented here in Appendix B). Figure |5| 
presents the calculated surface ionization probability of 
lithium for a residual air pressure of2.0xlCT 7 Torr, com- 
pared to the experimental results of Kawano et al. 0. 
The present microscopic results, obtained with the sur- 
face work function modeling, are seen to reproduce the 
experimental trends rather well, the abrupt change of the 
ionization probability from 100% down to a few % as the 
temperature is increased clearly appears to be connected 
with the change of surface work function or equivalently 
to the degree of oxidation of the rhenium surface. The 
two limits (high and low T) in the present case corre- 
spond to cither an almost clean rhenium surface and an 
oxidized rhenium surface. The temperature T c at which 
the abrupt change occurs depends on the residual gas 
pressure which directly influences the oxygen adsorption 
change with temperature. T c decreases when the residual 
pressure decreases. Typically, it changes by around 300 
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K for a residual air pressure change between 10 -9 and 
2 x 1CP 6 Torr. The experimental and theoretical vari- 
ations of the ionization probability are shifted one with 
respect to the other by around 40 K. This can be due 
to the approximations involved in the present modeling 
and/or to the accuracy of the temperature scale (see be- 
low and Appendix A). Finally, the asymptotic value of 
the ionisation probabibility at high temperature is under- 
estimated by our calculation and this proves that we use 
a too low value of the work function <& of clean unoxi- 
dized rhenium. Clearly, this value, <&o = 4.94 eV, taken 
from our fit of Kawano et al. data ( see Appendix 
B), is slightly lower than, for instance, the value obtained 
by Persky @, $ = 5.0 ± 0.1 eV. 



pearing from the result. One can notice that the freez- 
ing distance approximation ([)]) consists in applying the 
Saha-Langmuir law expression for the electronic levels 
alone at the point where the electronic levels decouple; 
the sum over the heavy particle energies (equation (^|)), 
which takes into account the energy changes between zf 
and infinity, transforms it into the usual Saha-Langmuir 
expression for the total energy of the system evaluated at 
infinity. The good agreement between the present micro- 
scopic model and Saha-Langmuir expression thus proves 
that dynamical effects are absent and that the specificity 
of the metal surface-projectile charge transfer process dis- 
appears. This comforts the validity of a thermal equilib- 
rium approach for the surface ionization process. 
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FIG. 5: Ionization probability of lithium on rhenium as 
a function of surface temperature with an air pressure of 
2 x 1CP 7 Torr: the full curve is the Saha-Langmuir law, us- 
ing our modeling of the work function (Appendix B); the 
empty triangles represent the theoretical calculation described 
in part 6, with the same value of the work function; the dots 
represent the data of Kawano et al. [Q] 

For surface ionization of sodium on rhenium surface, 
a similar agreement (not shown here) is found between 
the present modeling and the experimental results of 
Kawano et al. M. Figure || also presents the prediction 
of the Saha-Langmuir equation for the same conditions. 
The predictions of the present microscopic study are ex- 
tremely close to those of the Saha-Langmuir equation, 
typically within a couple of %. In fact, this can be un- 
derstood if we replace the value of from the numerical 
solution of equation (0) by its freezing distance approx- 
imation (281, i.e. if the final charge state P s + is taken 
equal to its equilibrium value at the freezing distance: 



P H 



2exp(^^; 



(9) 



Bringing equation (^|) into equation (||) then leads to the 
Saha-Langmuir equation ([!]), the value of E a (zp) disap- 



FIG. 6: Ionization probability of lithium on rhenium as 
a function of surface temperature with an air pressure of 
2 x 10" 6 Torr: the full curve is the our theoretical model- 
ing of the ionization process described in part 6, using the 
equations of Appendix B to describe the work function; the 
empty circles represent this same theoretical modeling calcu- 
lation, but using directly the measurements of work function 
[ fTi] ; the dots represent the data of Kawano et al. jj] and 
the diamonds our results (shown in figure ^|) multiplied by a 
factor 1/-E+ = 1.41 to get the true ionization probability. 

Figure || presents a comparison of experimental and 
theoretical results for the surface ionization probability 
for a residual air pressure of2.0xl0 -6 Torr. Two theoret- 
ical results are presented which were obtained either by 
using directly the surface work function extracted from 
electron emission experiments [|lT| or by using its mod- 
eling (Appendix B). The two results are close one to the 
other, confirming the efficiency of the modeling of the 
work function change. The two experimental sets are the 
results from reference and the present results. The 
latter have been multiplied by 1.41 = 1/P+ to transform 
our detection probability result D + into an ionization 
probability P + . Although the general behaviour of P+ 
as a function of T is the same in the three sets, they ap- 
pear shifted one with respect to the other. The abrupt 
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change in P + as a function of T from ref. M is more rapid 
than in the present study. In both cases, there is an up- 
ward temperature shift when going from the model to the 
experiments (about 70 K for the data of ref. j?J and 110 
K for our results). These differences are tentatively at- 
tributed to the accuracy of the temperature scales, both 
in the experimental results on surface ionization (esti- 
mated to 5 % in the present study) and in the modeling 
via the use of the electron emission experimental results 
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Appendix A: Hot wire temperature 



where / is measured in ampere. The numerical factor 
1124 corresponds to the dimensions (a = 30 pm and 
b = 760 pm) of the ribbon used in Toulouse and, for 
other ribbon geometries, this factor is easily scaled, using 
equation ([To|). All temperatures appearing in this paper 
were deduced from this equation (even up to 2250 K), 
without recalling our estimated 5% error bar. 

We have also used an optical pyrometer to measure 
the temperature: the sensitivity is close ±5 Kelvin, but 
we have no way to test its calibration. Moreover, such 
pyrometers are calibrated to measure the temperature of 
a blackbody radiation and the readings T r must be cor- 
rected to take into account rhenium spectral emittance 
Q near 655 nm C655 — 0.40 and window transmission 0. 
The corrected temperature T c is given by: 

1 c 1 r /If655 

This correction is substantial (about 200 K near 2000 
K). With our pyrometer, the readings T r are lower than 
the values of T deduced from equation (|l^) by about 5% 
while the corrected values T c are higher by roughly the 
same amount. 



For a long ribbon of thickness a and width 6, the equi- 
librium temperature T results from the equilibrium be- 
tween the input power due to Joule effect and radiative 
losses: 



p(T)P 



e(7>T 4 [2(a + &)] 



(10) 



where p{T) is the electrical resistivity, e(T) the to- 
tal emittance, both temperature dependent and a the 
Stefan-Boltzmann constant. It is a very good approxi- 
mation to forget thermal conduction along the wire, as 
long as one does not want to describe the temperature 
distribution near the ribbon ends. Rhenium resistivity in 
n ■ m JI§ is well fitted by: 

p(T) = 26.0 x 10~ 8 x (1 + 1.27 x 10" 3 x T) (11) 

for 1200 < T < 2000 K. For the total emittance of rhe- 
nium, we have also fitted the three data sets collected in 
ref. Q and covering the range 300 < T < 3000 K by a 
linear function of T: 



e(T) = 0.0852 x (1 + 1.15 x 10~ 3 x T) (12) 

all the data points being within ±20% of this fit: this un- 
certainty produces the dominant temperature error bar 
equal to 5%. It appears from equations ( [Ti"|]l2] ) that the 
ratio p(T)/e(T) is a very slow function of the temper- 
ature T. This explains why the simple law T cx \fl is 
an excellent approximation as noted by ref. ]l2[ ]. In the 
range 1200 - 2000 K, we get: 



Appendix B: Rhenium work function $ as a function 
of oxidation and temperature 

The strong influence of oxidation of rhenium surface on 
the ionization probability of lithium was first observed by 
Persky § in 1968. This study was continued by Kawano 
and coworkers jj], ||, (To), who used the Saha-Langmuir 
law to deduce the work function from the observed ion- 
ization probability. In 1999, Kawano et al. [|ll] also 
measured the work function of oxidized rhenium from 
the emitted electron current, using the Richardson law 
( J = AT 2 exp (— <fr/k B T)) to relate the current density 
J to the work function $. The two values of the work 
function differ noticeably, the work function for positive 
ion emission deduced from the Saha-Langmuir law being 
substantially larger than the work function for electron 
emission deduced from the Richardson law. 

We first discuss the work function for electron emis- 
sion because Kawano and coworkers |ll[] have collected a 
large data set as a function of oxygen pressure and rhe- 
nium temperature. We have developed a simple model 
which fits these results very satisfactorily. The oxygen 
surface coverage is assumed to be described by the fol- 
lowing chemical equilibrium: 

2 rhenium sites + O2 2 oxidized rhenium sites (15) 

Let a R e and <jR e o be the density of unoxidized and ox- 
idized sites respectively, and p(C?2) the partial pressure 
of molecular oxygen. The equation resulting from this 
chemical equilibrium is: 



T = 1124 x V7 



(13) 



VReO 



v{Qi) 

Pc(T) 



(16) 
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where p c (T) is the equilibrium constant of the reaction, 
written so as to have the dimension of a pressure. From 
this equation, we can deduce the fraction / of the oxi- 
dized sites: 

, = <?ReO = y/p{Q 2 )/Pc{T) 

ORe+OReO 1 + y/p(0 2 )/ Pc (T) 

We would get a different equation if we consider that 
each rhenium site accepts two oxygen atoms. We have 
also tried to fit the data with this modified equation. As 
the resulting fit is considerably less good, we consider 
that this second hypothesis is not correct. The second 
assumption of the model is that the work function $ 
increases linearly with the coverage / from the pure metal 
value $o to the completely oxidized value noted <£>o + A: 

$ = $o + /A (18) 




1 — i — i 1 — i — i — i 1 — i — i — i 

1x10" 7 1x10 s 1x10 5 1x10" 



Air pressure (Torr) 

FIG. 7: Work function of rhenium as a function of air pres- 
sure in Torr (logarithmic scale) for various temperatures: di- 
amonds for T — 1406 K, stars for T = 1527 K, triangles for 
T = 1685 K, dots for T = 1828 K and squares for T = 2015 
K. The experimental results are taken from figure 4 of the 
paper by Kawano et al. p3] , while the curves represent our 
best fit to these data, using equations ( fnl ) and (|lq). 

This assumption is surely oversimplified but we do not 
see how to refine easily this model as it already repre- 
sents quite well the results presented of H. Kawano et 
al. |ll| in the temperature range 1400 — 2000 K. Our 
fit to these results is presented in figure 0. Below 1400 
K, the observed variations of the work function deviate 
strongly from the trend observed at higher temperatures 
and we have not tried to fit this data. From the fits, we 
extract the following values: $ = 4.94 eV, A = 1.125 
eV and for each temperature a value of the equilibrium 
constant p c (T). The data of figure 5 of reference Q 
are well fitted too, but with a slightly lower A value, 
A = 1.00 eV. From these fits, we have deduced 8 values 
of p c (T). A good test of the consistency of our model is 



that, as expected from thermodynamics, these values are 
well represented by the following formula: 

log w (Pc(T))=A+^ (19) 

If the pressures are expressed in Torr, we get A = 9.186 
and B = —25072 K. The quantity B corresponds to an 
energy B = —2.16 eV, which represents the energy differ- 
ence between the two sides of the chemical equilibrium 
described by equation (|l^). This modeling provides a 
very efficient way of interpolating between the measured 
work function values and we used it in our treatment of 
the surface ionization to define the surface work function 
as a function of the operating conditions. 

The work function of rhenium can also be deduced 
from the measurements of the ionization probability of 
sodium or lithium (?|, ||, [l0| . The analysis is somewhat 
more complex because the experiments were done with 
intense alkali halide beams and the dissociation equilib- 
rium of the molecule on the surface must be taken into 
account. The Saha-Langmuir law is used to extract the 
work function <E> + (thus noted as it differs from the elec- 
tron emission value) from the data as a function of the 
surface temperature and oxygen pressure. However, two 
aspects of the analysis of Kawano and coworkers deserve 
further discussion: 

• the ionization potential / of the alkali atom is taken 
as a function of the wire temperature (this function 
appears to be I(T) = I + 2.5fcsT). This is a sub- 
stantial effect which is not commonly considered 
(see for instance ref. [ fL3] , in which very accurate 
measurements of the platinum and iridium work 
function are made as a function of the tempera- 
ture). Adding this unusual T dependence to the 
atomic ionization energy significantly contributes 
to the observed difference between the surface work 
function extracted from electron emission and the 
effective work function extracted from surface ion- 
ization data. 

• it is surely very difficult to estimate the ionization 
probability to better than a few percent, in partic- 
ular because the flux incident on the wire is calcu- 
lated from vapor pressure data which are not very 
accurate (as discussed in the present paper in the 
lithium case). It is therefore very difficult to under- 
stand how the work function $ + can reach a value 
as large as 6.7 eV (see figure 3b of ref. || where 
the experiment was done with NaCl): as soon as 
($ + — I) > 5.3/cbT, the Saha-Langmuir law pre- 
dicts P + > 0.99, almost impossible to distinguish 
experimentally from P + = 1. For temperatures 
close to 1500 K, the maximum value of <f> + , which 
can be reliably deduced from Saha-Langmuir law, is 
not larger than 6.1 eV when working with lithium 
(and only 5.8 eV when working with sodium for 
which / = 5.139 eV). 
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• finally, one can stress that most of the experimen- 
tal results of ref. j?J and || on the surface ioniza- 
tion probability of lithium and sodium can be re- 
produced by Saha-Langmuir law using the surface 
work function extracted from electron emission ]Tl| ] 
and the ionization potential of the free atom. There 
does not seem to be any need to introduce an ef- 
fective work function for surface ionization, which 



could be, at best, only a parameterization of the 
experimental results. 

Therefore, we think that the experimental results on 
P+, obtained by Kawano and co-workers are very inter- 
esting for the characterisation of a surface ionization de- 
tector, but the values of <£>+ extracted from these should 
be considered with caution. 
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